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Abstract. It is wcll known that, beginning in 2000, the behavior of the thermal 
correction to the Casimir force between real metals has been hotly debated. As was 
shown by several research groups, the Lifshitz theory, which provides the theoretical 
foundation for the calculation of both the van der Waals and Casimir forces, leads to 
different results depending on the model of metal conductivity used. To resolve these 
controvèrsies, the theoretical considerations based on the principies of thermodynamics 
and new experimental tests were invoked. We analyze the present status of the 
problem (in particular, the advantages and disadvantages of the approaches based 
on the surface impedance and on the Drude model dielectric function) using rigorous 
analytical calculations of the entropy of a fluctuating field. We also discuss the results 
of a new precise experiment on the determination of the Casimir pressure between two 
parallel plates by means of a micromechanical torsional oscillator. 



PACS numbers: 12.20.Fv, 12.20.Ds, ll.10.Wx, 05.70.-a 

1. Introduction: Lifshitz formula provides the theoretical background for 
the thermal Casimir force 

Since the beginning of 2000, the dependence of the Casimir force acting between real 
metals on separation and temperature has been hotly debated. It is well known that 
van der Waals and Casimir forces act between closely-spaced macroscopic bodies. Both 
forces are of the same physical nature and are caused by long wavelength electromagnètic 
fluctuations £Q . The name "van der Waals" is usually used in a nonrelativistic case when 
the force and other physical quantities depend on the Planck constant h, but do not 
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depend on the velocity of light c. In the relativistic case, which corresponds to larger 
separations, the name "Casimir" is commonly used, since he was the first to predict the 
force acting in vacuum between two paral·lel uncharged plates made of ideal metals 
During the last few years dispersion forces (which is the genèric name for both kinds 
of forces caused by fluctuations) were actively investigated, not only theoretically but 
also experiment ally (see the monographs [3-5] and a review (§]). This was motivated 
by many prospective applications of dispersion phenomena in fundamental physics and 
nanotechnology. In particular, the problem of an accurate theoretical description of 
dispersion forces between real bodies (i.e., taking into account the realistic conductivity 
properties, surface roughness and nonzero temperature) took on great significance. This 
paper is devoted to the outstanding problems arising in the case of plates made of real 
metals and having nonzero temperature. 

The free energy of the dispersion interaction between two thick parallel plates 
(semispaces) described by a dielectric permittivity e(u), with a gap of width z between 
them in thermal equilibrium at temperature T, was derived by Lifshitz 7J. It can be 
represented in the form 

= (i - 5*.) jf * A (1) 

X {lnfl-rjtó.tije- 2 »"] + ln [l - rifó, fcje" 2 ""] } . 

Here ks is the Boltzmann constant, = 2irkBTl/h are the Matsubara freqüències, 
qi = {k\ + £ 2 /c 2 ) 1//2 , k± is the projection of the wave vector in the plane of the plates, 
and riu (£,i,k±) are the reflection coefficients for two independent polarizations of the 
electromagnètic field. The explicit form of these coefficients in the case of real metals is 
closely connected to the controvèrsies which are the subject of this paper. 

2. Two main approaches to the presentation of reflection coefficients: the 
dielectric permittivity and the surface impedance 

In the framework of the first approach, the reflection coefficients are expressed in terms 
of the dielectric permittivity of a plate material. In application to Eq. (pQ) this approach 
goes back to E. M. Lifshitz who used it to describe the Casimir and van der Waals forces 
between dielectrics. In this case 

where k\ = {k\ + Si^f/c 2 ) 1 ^ 2 and E\ = e(iÇi). Within the Lifshitz approach, the 
fluctuating electromagnètic field is considered both in the gap and in the interior of 
the dielectric semispaces with the continuity boundary conditions on their surfaces 

E\ t = E 2t , B u = B 2 t, D ín = D 2n , B ln = B 2n (3) 

(note that we do not consider ferromagnets and set B=H). 

Of prime importance is the contribution of the free energy in Eq. (JIJ with l = 0. 
It is well known that this term determines the total value of T at high temperatures 
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(large separations) because all contributions with l > 1 are exponentially small in this 
limiting case j3Ej- Eq. (JTJ) with reflection coefficients (J2J) correctly describes the case 
of plates made of ideal metals if one uses the so-called Schwinger prescription JHJ IHJ, 
i.e., take limit e — > oo first and set Z = afterwards. In doing so for ideal metals the 
equations 

rj(0,k±)=rï(0,k±) = l (4) 

are obtained which coincide with those obtained from thermal quantum field theory 
with Dirichlet boundary conditions [H] independent of the Lifshitz formula. 

In Ref. [9j an attempt was undertaken to describe the dispersion interaction between 
real metals by the Lifshitz formula using the Drude model to characterize their dielectric 
properties at low freqüències: 



where u p is the plasma frequency and 7(T) is the relaxation parameter. Substituting 
Eq. (JHJ) into Eq. (J2J) one obtains 

rf (0, fcx) = 1, rí(0,fcj_) = 0. (6) 

Eq. (JHJ) is preserved also in the limit of ideal metals, and is thus in contradiction with 
Eq. (J3J). If, however, one characterizes a metal by the plasma model, 

em = l + f, (7) 
and extrapolates it to low freqüències, the reflection coefficients become: 



lc 2 k 2 , + u 2 — ck± 

ri(0,k ± ) = l, r ±(0> k ±)= \ 21 t \ r- (8) 

!C Z K\ + Lúí + CK\ 



Eq. (JHJ) agrees with the case of ideal metals (HJ) because r^(0, k±) — > 1 when uj p — > oo. 

The approach to the theoretical description of the thermal Casimir force between 
real metals based on Eq. (jüj) was supported and developed in Refs. [10-13], while that 
following Eqs. (J7|) and (JHJ) was proposed in Refs. ^2J EI] • By convention we refer to them 
below as "the Drude model approach" and "the plasma model approach" , respectively. 

The Drude model approach predicts relatively large thermal corrections to the 
Casimir force between real metals at short separations of about 500 times greater than 
between ideal metals. At large separations or high temperatures the Drude model 
approach predicts 1/2 the magnitude of the Casimir free energy and force than between 
ideal metals. (The same prediction was obtained using the model of a conducting wall 
by the classical Coulomb fluid fBJ or "non-relativistic quantum electrodynamics" |17j.) 
This prediction is in contradiction with the classical limit based on Kirchhoff's law [18J. 
As to the plasma model approach, its predictions are in qualitative agreement with the 
case of ideal metals. Because of this, the plasma model approach is called "traditional" 
and the Drude model approach is called "alternative" . The apparent advantage of the 
Drude model approach is that Eq. (JSJ) presents the same behavior of e ~ u)~ x at low 
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freqüències, as given by the Maxwell equations, whereas the behavior e ~ u~ 2 in Eq. (J7|) 
is characteristic for the region of infrared òptics, raising a question on the correctness 
of its extrapolation to low freqüències. On the other hand, Eq. (JHJ) is in agreement with 
the limiting case of an ideal metal, whereas Eq. (jüj) is not. Bearing in mind that the 
description of real metals using e depending only on u is not universal, it is worthwhile 
to use alternative physical quantities to characterize the reflection properties of metal 
surfaces. 

In the framework of the second approach, which goes back to V. L. Ginzburg and 
his collaborators (see, e.g., Ref. JH]), the reflection coefficients in the Lifshitz formula 
are expressed in terms of the Leontovich impedance Z(u). The Leontovich impedance 
boundary condition 

E t = Z(lu) [H t x n) (9) 

(where n is an infernal unit vector directed into the medium and normal to the 
boundary) is approximate and permits one to find the electromagnètic field in the 
gap without considering it in the interior of the plates. In terms of the Leontovich 
impedance, the reflection coefficients in Eq. are given by 

ni=>fte,*J = ^!!, rx S rf(6,tJ = ^^, (10) 
11 cqi + Zfa Íi + cqiZi 

where Z\ = Zii^i). The Leontovich impedance is applicable even in some cases when 
e(uj) loses its meaning (e.g., in the frequency region of the anomalous skin effect). 
If, however, both quantities are in their domain of application, they are related by 
Z(üü) = 1/ 'J ' s{uj) [2H]. For the impedance functions of the normal and anomalous skin 
effect, the Ginzburg reflection coefficients at zero frequency satisfy Eq. (JU) as do the 
Lifshitz ones for ideal metals. In the region of infrared òptics it follows that 

rf(0,^) = l, r G ± (0,k ± ) = ^^. (11) 

The impedance approach to the thermal Casimir force was developed in Refs. [21-24]. 
It predicts small thermal corrections at short separations in agreement with the case 
of ideal metals and leads to the same high-temperature (large separation) asymptotic 
results which hold for the plasma model approach or for ideal metals in accordance with 
the classical limit. Thus, the impedance approach can also be called "traditional" . 

Many computations of the thermal Casimir force have been performed in the 
framework of the Drude model approach [10-13], plasma model approach f3J ES] and 
impedance approach [21-24], and many controversial statements in favor and against 
each approach have appeared in the literature. These arguments were most recently 
presented in Refs. |2*ïï| I2fij. Before discussing them in more detail, we emphasize that 
the Drude model and the impedance approaches are qualitatively different only at 
zero Matsubara frequency [compare Eqs. (jüj) and (JHJ)]. As to the contributions from 
Matsubara freqüències with l > 1, both approaches find these using tabulated optical 
data extrapolated to low freqüències by the imaginary part of the Drude dielectric 
function with practically coinciding results. The plasma model approach applied at all 
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freqüències disregards the internal photoelectric effect (interband transitions) and other 
processes taken into account in the optical data. Because of this, it is not as accurate 
as the impedance approach at short separations of about 150-250 nm. 



3. Controvèrsies between the two approaches 



3.1. Is there a violation of thermodynamics for the Lifshítz formula combined with the 
Drude model? 



As was proven analytically in Ref. j23|, the Drude model approach leads to a violation 
of the third law of thermodynamics (the Nernst heat theorem) in the case of metallic 
perfect lattices with no defects or impurities. For such lattices the relaxation parameter 
7(T) — ► when T — > in accordance with the Bloch-Grüneisen law (this property is 
preserved even when the effects of electron-electron col·lisions are included), and the 
entropy of a fructuating field takes the form 

2 



S(z,0) 



dT{z, T) 



k 



B 



dT 



T=0 



IQtxz 2 



ydy ln 



1 



y 



< o, 



(12) 



where uj c = c/(2z) is the characteristic frequency of the Casimir effect. 

To avoid this conclusion, Refs. ^21 EU] apply the Drude model approach to 
metallic lattices with defects and impurities possessing some residual relaxation 7(0) ^ 
and obtain S(z, 0) = 0. This, however, does not solve the problem of the 
thermodynamic inconsistency of the Drude model approach, because perfect lattices 
have a nondegenerate dynamical state of lowest energy and, thus, according to quantum 
statistical physics, the entropy at zero temperature must be equal to zero in this case 
[a property violated by the Drude model approach according to Eq. (fT2"|) ]. Recently the 
authors of Ref. [26J recognized that the Drude model approach violates thermodynamics 
"for perfect metals of infmite extension" having no relaxation (the so-called "modified 
ideal metal" which is obtained from the real metals in the Drude model approach when 
e — ► 00 but they deny a violation for metals which include relaxation. This 

denial is based on a misunderstanding. As was already mentioned above, the Drude 
model approach violates the Nernst heat theorem for metallic perfect lattices with no 
defects or impurities [see Ref. j2lj and Eq. (fï2"j) ]. Such lattices have a nonzero resistance 
and relaxation at any nonzero temperature, which go to zero only when temperature 
vanishes. Thus, they are not "perfect metals" according to Ref. [26J. 

In the light of the above, the conclusion of Refs. ^21 EH] that the Drude model 
approach is in agreement with thermodynamics cannot be supported. In contrast, 
the plasma model and the impedance approaches are in complete agreement with 
thermodynamics fSJ 121] . 
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3.2. Do the "exact" impedances lead to the same result as the Drude model or as the 
Leontovich impedance? 

As was mentioned in Sec. 2, the impedance boundary condition (jHJ) is an approximate 
one. In particular, it works well if the electromagnètic waves penetrating into metal 
interior propagate almost perpendicular to its surface. (This is in fact the long 
as \e(uj)\ 1.) Some authors (see, e.g., Refs. O 123) introduce so-called "exact" 
impedances 

(2t,2 \ / 2r,2 \ — 1/2 

, Z ± (u,k ± ) = Z(u)íl-^j , (13) 

which depend on both polarization and transverse momentum, and related impedance 
boundary conditions 

E t = Z±(u, k±) [H t x n] , Z|| (oo, k ± )H t = [n x E t ] , (14) 

where Z(u) is the Leontovich impedance. Actually, the impedances (fTïïj) are not exact. 
Eqs. (|13|) and (|14j) follow from the Maxwell equations only if e(u) has definite meaning, 
so that the spatial dispersion is neglected and the relation D{r,oj) = e{uj)E{r,uj) is vàlid. 
If, however, the homogeneity of space is violated by the presence of boundaries, and if 
spatial dispersion is taken into account, Eqs. (|T3*j) and ()14j) are in general inapplicable 
(see Sec. 5). 

In terms of the impedances (|13|) . the reflection coefficients take the form 

cqi + Zi^i 4; + cqiZ ií± 

where Z t j = Z\\(iÇi, k±) and Z^± = Z±(iÇi, k±). It is easy to check that if one substitutes 
Eq. (fTïïjl into Eq. (JT3J), considers £i and k± as independent variables, and takes into 
account that Z = 1/ y/e, the Lifshitz reflection coefficients (j2J are recovered. This 
allowed the authors of Refs. jTTl[T2||2n] to insist that the results (HJ) and (JTTJ), obtained 
in Refs. P2^l2*4^ within the impedance approach [in contrast to Eq. (jHJ) vàlid in the Drude 
model approach], are due to disregarding the transverse momentum dependence by the 
Leontovich impedance. 

Bearing in mind that Eq. (jHJ) returns us to the contradiction with thermodynamics, 
it is necessary to analyze this point more thoroughly. When we are considering real 
photons incident on a plane boundary of a single semispace, the components of a wave 
vector are constrained by the mass-shell equation \k\ 2 = k\ + k\ = uj 2 /c 2 (k 3 is the 
wave vector component perpendicular to the boundary). Using this equation, the angle 
of incidence is expressed by sin^o = k±/\k\ = ck^/u and the impedances (fTïïjl can be 
identically rewritten in the form 



(w,k\) = Z(üj) 



sin 2 $r 



e{u>) 



1/2 



Z±(uj, k±) = Z(uj) 



sin 2 -d 
e(uj) 



-1/2 



(16) 



The term sin 2 $ /e(uj) can be neglected for all freqüències which are at least several 
times smaller than the plasma frequency because in this region |e| ^> 1. Furthermore, 
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in the limit of zero frequency sin 2, ü /e — > due to |e| — > oo, and as a result for real 
photons the Leontovich impedance coincides with the impedances Zn± precisely. The 
mass-shell equation ensures the fulfilment of the Snell's laws. Thus, if the reflection 
properties for virtual photons are postulated to be the same as for real photons, the 
Leontovich impedance and boundary condition (jHJ) are recovered. 

These considerat ions can be generalized for the case of two semispaces with the 
boundary conditions (JHJ). 111 this case the Maxwell equations lead to the following 
dispersion relations for the determination of photon eigen-frequencies uj [22J, 

A|| (uj, k ± ) = e~ qz ( sinh qz - cos h qz } = 0, (17) 



where 




2ÍK± \ 

A±(uj, k±) = e~ qz ( 1 — Kj_ ) [ sinh qz + - ^ cosh qz \ = 0, 



K 



2 ,2 2/2 uZ\\(uj,kA_) cqZ ± (uj,k±) 

q =k ± -u/c, r?|| = ü , k±= . (18) 

cq uj 

It is easily seen from Eq. (JT7|) that in the limit uj — > the transverse momentum must 
also go to zero, k± — > 0, in such a way that c/c^ ~ c<j. This means that the quantity 
c 2 kj_/ {uj 2 e) in Eq. (JT3j) goes to zero in the limit of zero frequency, and both impedances 
coincide precisely with the Leontovich impedance. At all nonzero freqüències there are 
only negligible differences between the impedances and the Leontovich impedance. 
Thus, in the case of two semispaces, the dispersion relations (fTTj) play the same role as 
the mass-shell equation for free photons. 

Using the above property of the solutions of Eq. ([T7)l. this equation can be 
approximately rearranged to the form following from the boundary conditions (jOJ) , 

A|| (u, k±) = e- qz (l - r/ 2 ) ( sinh qz - ] — cosh qzj = 0, (19) 




A_i_(cj, k±) = e qz (l — k 2 J ^sinhgz + - — ■ — - cosh qz^j = 0, 

where rj = ujZ/ (cq), k = cqZ/uj are now expressed in terms of the Leontovich impedance. 
Note that the approximation was made only at uj ^ 0, whereas in the limit uj — > 0, which 
is of special importance for our problem, Eqs. (fPTj) and ()19j) remain equivalent. 

It is now easy to proceed with the derivation of the Lifshitz formula by summing 
the free energies of all oscillators whose freqüències are determined by equations í)19|l. 
applying the argument theorem and calculating the residues at imaginary Matsubara 
freqüències [22] • This derivation results in Eq. (JI} with the reflection coeficients in 
Eq. pOjl . Note that in our derivation the transverse momentum dependence is not 
"completely disregarded" (as is claimed in Refs. fU^H^E])- The transverse momentum 
k± enters through the quantity q both in the dispersion equations (fTïïjl and in the 
reflection coefhcients fil 0(1 . We emphasize that the dispersion relations f|19j) with the 
Leontovich impedance have both photonic and plasmonic solutions. It is not admissible, 
however, to preserve k± ^ in the impedance function Z±(uj, k±), and to consider uj — > 
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because this violates the dispersion relations (|17|). The Matsubara freqüències and 
the transverse momentum become independent only after application of the argument 
theorem. In view of the above, the conclusion of Ref. ^2] that the Drude model and 
impedance approaches lead to the same reflection coefficients at zero frequency is not 
warranted. 

3.3. How to extrapolate real data to zero frequency? 

As was already discussed in Sec. 2, the tabulated optical data for real materials 
do not contain the zero frequency contribution which should be found from some 
theoretical considerations. There are two main methods proposed in literature for 
solving this problem. According to one method (see, e.g., Refs. fT2". 26J) it is necessary 
to extrapolate to zero the actual reflection properties of plate materials at very low, 
quasistatic freqüències. According to another method, the zero-frequency limit should 
not be understood literally because there are no static-field fluctuations. More likely 
it should be understood as a mathematical limit to zero from the region around the 
characteristic frequency uj c giving the major contribution to the thermal Casimir force 
[21-24]. Although in Ref. [12] such an extrapolation is considered to be "ad hoc", in 
Ref. j2Hj it was demonstrated that the method of Refs. [121 EH], applied in the framework 
of the impedance approach, results in a violation of the Nernst heat theorem. 

Recently the method of Refs. [21-24] received support when applied to two dielectric 
semispaces, i.e., to the case which was previously free of any controvèrsies. As was 
analytically proved in Ref. [2H1 fo r similar (and in Ref. [29J for dissimilar) dielectrics, 
the Casimir entropy of a fluctuating field goes to zero when temperature vanishes if the 
semispace materials are described by a frequency dependent e with some finite static 
vàlues e(0). If, however, the dc conductivity of dielectrics is taken into account, this 
results in a violation of the Nernst heat theorem. When it is recalled that dielectrics 
really possess some nonzero conductivity at a constant current (although it is many 
orders of magnitudes smaller than that for metals), it becomes clear that real material 
properties at very low freqüències are in fact irrelevant to the fluctuating phenomena 
and must not be included into the model of the dielectric response. 

4. Different approaches to the thermal Casimir force and experiment 

4-1. Measurement of the Casimir force using a torsion pendulum 

As was mentioned in Sec. 2, the traditional approaches to the thermal Casimir force 
predict small thermal corrections at short separations in qualitative agreement with the 
case of ideal metals. The sensitivity of all already performed experiments is not sufficient 
to measure such corrections, much less to discriminate among different traditional 
approaches. Some of these experiments are, however, sufficiently sensitive to detect 
the alternative thermal corrections as predicted by the Drude model approach if they 
exist. 
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According to the results of the first modern measurement of the Casimir force 
between a Au coated plate and large spherical Iens jSÜ] , at T = 300 K with z = 1 /xm, 
a net deviation between the Casimir forces computed for ideal metals and using the 
Drude model approach is 25%. Of this deviation, 19% is due to the large alternative 
thermal correction. In Ref. [HH) the experimental uncertainty at z — 1 /im was about 5- 
10%. The predicted thermal effect was not, however, observed. According to Ref. jÏÏT], 
the prediction of the Drude model approach "deviates significantly from experimental 
results". At the same time, according to Ref. [H2j, the experimental data are consistent 
with the impedance approach. 



4-2. Static experiment using a microelectromechanical torsional oscillator 

In this experiment the Casimir force between a plate of a torsional oscillator and a 
sphere was measured [33] . The sphere was coated by Au and the plate by Cu. The 
force vàlues were measured with an absolute error of 0.3 pN for separations in the 
range 0.19-1.2 /im. Ref. |T2] claims that there is "reasonably good agreement" between 
the Drude model approach and the experimental results of the static experiment at 
z = 200 nm. In Ref. [22] it was demonstrated, however, that this conclusion is based 
on a misunderstanding of the data. In fact, the experimental data are in contradiction 
with the Drude model approach, but are consistent with the traditional approaches 
to the thermal Casimir force. For example, at about z = 200 nm the mean value 
of F th — F cxp (where the theoretical Casimir force F th is computed using one of the 
traditional approaches) is practically zero. If the Drude model approach is used, this 
difference is approximately 2.6 pN (i.e., it is almost 9 times larger than the experimental 
absolute error), thus, upsetting the agreement between experiment and theory. 

4-3. Dynamic experiment using a microelectromechanical torsional oscillator 

In the second, dynamic, experiment of Ref. jHSj the equivalent Casimir pressure 
between two parallel plates (one coated by Au and another by Cu) was determined 
dynamically in the separation range from 260 to 1200 nm with an absolute error of about 
0.6 mPa. In the comparison between the experimental data and different theoretical 
approaches, a number of possible corrections to the theoretical Casimir pressures were 
taken into account. These include those due to surface roughness, the use of the 
proximity force theorem, the finite sizes of the plates, sample-dependent variations of the 
tabulated optical data for the complex index of refraction, and contribution of nonlocal 
effects. The differences between the theoretical and experimental Casimir pressures, 
P th (z) - P cxp (z) , deviate from zero in the range from -0.8 to 0.5 mPa if the traditional 
theoretical approaches to the thermal Casimir force are used (3S]- This demonstrates 
that the traditional approaches are consistent with the experimental data. If, however, 
the theoretical Casimir pressures P th (z) are computed using the Drude model approach, 
the differences |^P th (z) — P cxp (z) deviate significantly from zero within the separation 
region from 260 to 700 nm. At the shortest separation this deviation reaches 5.5 mPa 
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(i.e., it is more than 9 times larger than the experimental absolute error). In Ref. [33] 
the conclusion was drawn that the Drude model approach is ruled out experimentally. 

4-4- Improved dynamic experiment using a microelectromechanical torsional oscillator 



In Refs. [34) 135] the results of a new, improved, dynamic experiment were reported on 
the determination of the Casimir pressure between the two plates both coated by Au in 
the separation region from 160 to 700 nm. The main improvements were a significant 
suppression of the surface roughness, a decrease by a factor 1.7 in the error in the 
measurement of the absolute separations, and the use of rigorous statistical procedures 
in data processing and in the comparison of experiment and theory. Importantly, the 
contribution of the surface roughness to the Casimir pressure was reduced to less than 
0.65% even at the shortest separation. This is the first experiment where the total 
relative error of the Casimir pressure measurements found at 95% confidence varies 
between 0.54 and 0.59% in a wide separation region from 170 to 300 nm, and is as small 
as 2.5% at z — 500 nm. 

In addition to the above-mentioned corrections in the theoretical Casimir pressure, 
the contribution of patch potentials and errors in the calculated pressures, arising 
from the uncertainties in the experimental distances |36 j were taken into account. 
Finally, the total theoretical error, which was determined at 95% confidence, decreases 
from 1.65% at z — 160 nm to 0.9% at z = 750 nm. Using the total experimental 
and theoretical errors, determined independent ly, the absolute error of the pressure 
differences P th (z) — P cxp (z) was found at 95% confidence. It is practically the same 
in different theoretical approaches. The sòlid lines in Fig. l(a,b) show the confidence 
interval versus separation which is determined in terms of this error. In Fig. l(a) the 
differences P th (z) — P exp (z) for the traditional theories are plotted. In this case more 
than 95% of all individual points (not only the mean vàlues of pressure differences at 
each separation as is demanded by the theorems of mathematical statistics) belong to 
the confidence interval. This means that our error analysis is in fact too conservative and 
the errors are significantly overestimated. From Fig. l(a) it follows that the traditional 
approaches to the thermal Casimir force are consistent with data. In Fig. l(b) the 
theoretical pressures P th (z) are computed using the Drude model approach. It is seen 
that this approach is excluded experimentally within the separation region from 170 to 
700 nm at 95% confidence. As is shown in Ref. [31], in the separation region from 300 
to 500 nm the Drude model approach is excluded by experiment at even higher 99% 
confidence. 

In the first version of preprint |26j . it was claimed that the computations of Ref. [31] 
in the framework of the Drude model approach deviate somewhat from the original 
computations of Ref. [33- Later it was shown [23] that these deviations were caused by 
numerical discrepàncies in the data used in preprint [26 , which was recognized by the 
authors in the second version of this preprint. 

Nevertheless, Ref. j2S] claims that "the experimental situation is still too indecisive 
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Figure 1. Differences between the theoretical Casimir pressures calculated using one 
of the traditional approaches (a) or the Drude model approach (b) and experimental 
ones versus separation. The 95% confidence interval lies between the sòlid lines. 

to draw definite conclusions" on the contradiction of the Drude model approach with 
experiment. To justify this conclusion, Ref. (2Ej cites Ref. [SE] where the influence of the 
experimental uncertainty in surface separations on the theoretical vàlues of the Casimir 
force was discussed. As was mentioned above, the full contribution from this effect is 
included in the error analysis of Ref. [33j , and thus cannot be used to cast doubts on the 
accuracy of the obtained results. Ref. also claims that "the roughness of the surfaces 
is much larger than the precision stated in the determination of the separation" . This, 
however, cannot compromise the precision achieved in Ref. [33] because the separations 
are measured between zero levels of surface roughness independently of the value of 
roughness amplitude. Furthermore, the contribution of roughness in the experiment of 
Ref. jH3] is so small, that the validity of all the conclusions is preserved even if one were 
to neglect the roughness in the theoretical analysis. As was demonstrated in Ref. |25j . no 
unaccounted fixed systematic error in surface separation is capable of bringing the Drude 
model approach into agreement with data within the whole range of measurements. One 
more claim of Ref. is that "accurate determination of a small difference between 
experimental vàlues at room temperature and purely theoretical vàlues at T = gives 
rise to further difficulties". This claim is evidently misleading because all theoretical 
vàlues of the Casimir pressure in Ref. [H3] were computed at T = 300 K, i.e., at the 
experimental temperature. To summarize, all of the arguments of Refs. ^21 EHj against 
the precision and accuracy of the experiment of Ref. [33] break down, and the conclusion 
that there is a contradiction between the Drude model approach and experiment seems 
unavoidable. 
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5. Spatial dispersion and the Casimir force 

Recent papers jHZl EH] claim that taking account of the spatial dispersion is crucial 
for the problem of the thermal Casimir force between real metals. They use the 
usual (continuity) boundary conditions (JHJ) and the usual Lifshitz formula (|TJ with the 
reflection coefficients expressed it terms of the spatially nonlocal dielectric permittivities 
Eí(uj, kx) or impedances Z»±(cü,k±) (the latter are presented in terms of nonlocal 
dielectric permittivities by means of some integral relations). Refs. [HZ1 EB] arrive 
at the conclusion that the contribution from r±(0,k±) is negligibly small, leading to 
approximately the same results for the thermal Casimir force as in the Drude model 
approach. In contrast to the Drude model approach, entropy was shown to be equal to 
zero at T = even for perfect crystal lattices [38J. 

As was already demonstrated in Sec. 4, the predictions of the Drude model approach 
are in contradiction with experiment. There are also serious theoretical objections 
against the formalisms used in Refs. [SUES]- It has been known that in the presence of 
spatial dispersion the continuity boundary conditions must be generalized for [SÜ] EO] 

47T 

E u = E 2t , B ln = B 2n , D 2n - D ln = Aïra, [n x (B 2 - B x \ = —i, (20) 
where 

li 2 li 2 dD 

, = -/div[„x[Dx„]]d/, i=^J lW dl (21) 

are the induced charge and current surface densities (assuming there are no external 
charges and currents). If spatial dispersion is absent, then a = 0, i— 0, and the usual 
continuity boundary conditions (jHJ) are recovered HO] • Then the use of boundary 
conditions (j3J) in the presence of spatial dispersion is unjustified. 

Refs. (HZ1 EH] disregard the fact that spatially nonlocal dielectric permittivities 
Eí(uj, kj_) can be rigorously introduced only for an unbounded medium which is uniform 
in space pïïH |4*U] . In the theory of the anomalous skin effect, such permittivities 
are sometimes used in the presence of boundaries (see, e.g., Ref. [41]). This is an 
approximation applicable to the investigation of some bulk effects, when physical 
phenomena caused by a layer adjacent to the boundary surface are neglected. (There 
is another approach to the theory of the anomalous skin effect [12] which only deals 
with the Leontovich impedance.) In any case, it is unlikely that the approximate 
phenomenological approach using nonlocal £í(u, k±) in the presence of boundaries would 
be applicable in the theory of Casimir force between metal surfaces, where the boundary 
effects are of prime importance. 

It has also been known |43 j that with inclusion of spatial dispersion the free energy 
of a fluctuating field T{z, T) takes a more general form than is given by the Lifshitz 
formula, i.e., 

F{z,T)=r(z,T)+AF(z,r), (22) 

where T{z ) T) is the Lifshitz contribution presented in Eq. í[T|). and the general 
expression for AjF(^, T) in terms of the thermal Green's functions of the electromagnètic 
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field can be found in Ref. [121 • The use of the Lifshitz free energy T{z ) T) for 
systems with spatial dispersion (as was done, e.g., in Ref. [221 where £j(o;,fcj_) of the 
hydrodynamic model were substituted into J 7 ) is not correct (23! • As was emphasized 
in Ref. (23j, "For most of condensed matter bodies this is inadmissible". It seems, 
however, that Refs. (33 EH! re P ea t this misuse of the Lifshitz formula. At the same 
time, it is amply evident that at T = 300 K in the region of infrared òptics (i.e., at 
the experimental separations) or in the region of the normal skin effect (z > 4 — 5 /im) 
spatial dispersion does not play any role and can be neglected. 

6. Conclusions 

From the above the following conclusions can be formulated: 

- The thermal Casimir force between real metals presents serious problems, and 
there is presently no fundamental theory which describes it; 

- The Drude model approach to the thermal Casimir force violates the Nernst 
heat theorem in the case of metallic perfect crystal lattices with no impurities and is 
excluded by experiment; 

- The thermal effect in the Casimir force has not yet been measured. The 
traditional approaches to its description based on the plasma model and on the 
Leontovich impedance are consistent with experiment. Computations using the Lifshitz 
formula at T = and tabulated optical data are also consistent with experiment; 

- The measurement of the small thermal effect, as predicted by the traditional 
approaches, would be of great interest and could help to discriminate between them. In 
this regard the proposed experiments [32 133 12S1 are very promising; 

- The desired fundamental theory of the thermal Casimir force between real 
materials should be presumably based on a more sophisticated quantization procedure 
incorporating spatial dispersion and inhomogeneity of space, and may go beyond the 
scope of the Lifshitz theory. 
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